Designing heterostructures with higher temperature superconductivity 
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We propose to increase the superconducting transition temperature T c of strongly correlated ma- 
terials by designing heterostructures which exhibit a high pairing energy as a result of magnetic 
fluctuations. More precisely, applying an effective theory of the doped Mott insulator, we envisage 
a bilayer Hubbard system where both layers exhibit intrinsic intralayer (intraband) d-wave super- 
conducting correlations. Introducing a finite asymmetry between the hole densities of the two layers 
such that one layer becomes slightly more underdoped and the other more overdoped, we evidence 
a visible enhancement of T c compared to the optimally doped isolated layer. Using the bonding and 
antibonding band basis, we show that the mechanism behind this enhancement of T c is the inter- 
band pairing correlation mediated by the hole asymmetry which strives to decrease the paramagnetic 
nodal contribution to the superfluid stiffness. For two identical layers, T c remains comparable to 
that of the isolated layer until moderate values of the interlayer single-particle tunneling term. These 
heterostructures shed new light on fundamental questions related to superconductivity. 

PACS numbers: 74.78. Fk, 74.20.-z, 74.72.-h 



I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity [T| , considerable efforts have been devoted to find- 
ing out how and why it works P-fTTj]. This puzzling 
phenomenon — electrical conduction without resistance 
at temperatures of up to ~ 130 K — occurs in com- 
plex "copper-oxide" materials (cuprates). After 1987, the 
term high-Tc superconductor was used interchangeably 
with cuprate superconductors until iron-based supercon- 
ductors were discovered [l2j, [l3| . 

The Hubbard model is a well-known model of in- 
teracting particles in a lattice, with only two terms in 
the Hamiltonian: a kinetic term allowing for tunneling 
( "hopping" ) of particles between sites of the lattice and a 
potential term consisting of an on-site interaction. There 
are many reasons to believe that the Hubbard model con- 
tains most (but maybe not all) of the ingredients nec- 
essary for understanding high-temperature superconduc- 
tivity 14j]. At zero hole doping, the single-band Hubbard 
model definitely captures the insulating behavior of the 
parent cuprate compounds. The origin of this insulating 
behavior was described many years ago by Nevill Mott 
as a correlation effect [l5j] and there is a suppression of 
the quasiparticle spectral weight 16]. In the Mott phase, 
electron spins form an antifcrromagnetic arrangement as 
a result of the virtual hopping of the antiparallel spins 
from one copper ion to the next — the parallel configura- 
tion being disallowed by the Pauli exclusion principle. It 
is relevant to observe that copper-oxide materials are gov- 
erned by a relatively large magnetic exchange J ~ 1300 
K which is much larger than the Debye energy of copper. 

Upon doping with holes the antiferromagnetism be- 
comes rapidly destroyed and above a certain level super- 
conductivity occurs with d x 2_ y 2 pairing symmetry. The 
d x 2_ y 2 wave nature of the order parameter has been con- 



clusively shown using phase sensitive experiments [17Hl| 
for example. The earliest experimental observation for d- 
wave symmetry is based on the linear decrease of the su- 
perfluid stiffness with temperature [2(|. An anisotropic 
gap with a d-wave order parameter has also been ob- 
served through photoemission studies [21!, 22]. One sce- 
nario to explain the d x 2_ y 2 wave nature of the supercon- 
ducting gap relies on spin fluctuations at the wavevector 
(jr, 7r) which makes the singlet channel attractive at large 
momentum transfer. This essentially stems from band- 
structure nesting effects in two dimensions close to half- 
filling [23]. A similar pairing occurs in ladder systems as 
a result of short-range valence bond correlations [2J, |25| ■ 
In 1986 it has also been suggested that backscattering 
from spin fluctuations might lead to the pairing seen in 
the Bechgaard salts [26| . The same year, three papers 
argued that spin fluctuations are resp onsible for d-wave 
pairing in heavy fermion systems [27H29| . 

How the electronic structure evolves with doping from 
a Mott insulator into a d-wave superconductor is a key 
issue in understanding the cuprate phase diagram. Over 
the years it has become clear that states in different parts 
of momentum space exhibit quite different doping depen- 
dences. The Fermi arcs [3(| EJ or pockets [32, 33] (near 
nodal states) retain their coherence as doping is reduced, 
while the antinodal (near the edge of the first Brillouin 
zone) states diminish in coherence, becoming completely 
incoherent at strong underdoping. The antinodes open 
a gap, "the pseudogap" , which appears well above the 
superconducting state. It is important to note that the 
relationship between the pseudogap and superconductiv- 
ity is still an open subject [3J, |35[ even though some 
efforts have been accom plis hed from the theoretical and 
numerical fronts [Til [36l445| ] . 

There are theoretical indications that the high T c in 
the cuprates may result from the large magnetic exchange 
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J @, H^-[13 • Designing a material that can increase T c 
certainly requires a better understanding of the mecha- 
nisms that reduce the superfluid stiffness with tempera- 
ture and with the proximity to the Mott insulating state. 
A number of recent theoretical [49l - l5l| and experimental 
[52I [53l | proposals have explored the possible benefit to 
combine quite metallic layers with layers of underdoped 
cupratc materials in heterostructured geometries. 

In this paper, we propose to increase T c using 
two strongly-correlated Hubbard layers, one layer be- 
ing slightly underdoped and the other rather overdoped; 
here, both layers are characterized by prominent d-wave 
correlations. Using an effective theory of the doped Mott 
insulator 0, H3] for both layers, we report an enhance- 
ment of the superconducting transition temperature com- 
pared to the optimally doped single layer. Another pos- 
sibility to increase T c relies on the presence of a very 
overdoped "free electron" like layer [49j (see Fig. 6). 

More precisely, using the bonding and antibonding 
band representation in the vicinity of optimal doping, 
we show that the low-energy BCS Hamiltonian exhibits 
dominant intraband d-wave pairing. Then, we justify 
how T c can be enhanced as a result of (interband) addi- 
tional superconducting fluctuations mediated by the hole 
asymmetry between the layers. Our results presented in 
Fig. 1 indeed reveal an enhancement of T c by ~ 20% 
for a slightly underdoped layer with hole density ^0.15 
and an overdoped layer with hole density ~ 0.25 — these 
include the possible charge redistribution when coupling 
the layers; optimal doping here refers to hole densities 
~ 0.17. Our findings may have applications to multilayer 
materials as well as heterostructures. The analysis per- 
formed in this paper uses a purely homogeneous model 
which does not include stripe or density wave structures 
53]. Other Hubbard bilayer systems have been stud- 
ied in different parametric regimes [HI, [HI- It is also 
worth mentioning that correlated bilayers exhibiting ei- 
ther heavy fermions I56l458j , composite fermions [591 ] or 
exciton condensates [60] are also attractive subjects. 

The realization of high-Tc superconductivity confined 
to nanometre-sized interfaces has been a long-standing 
goal because of potential applications [61 1. 

The outline of the paper is organized as follows. In 
Sec. II, we introduce the low-energy theory of the bi- 
layer system including the effect of Mott physics (large 
interactions) and we discuss the general methodology. In 
Sec. Ill, we address the situation of symmetric layers and 
show that there is no enhancement of T c ; nevertheless, 
we like to emphasize that the d-wave superconducting 
state is quite robust toward the proliferation of quasi- 
particles favored by the single-(quasi)particle tunneling 
term between the layers and therefore T c remains almost 
constant until moderate values of the interlayer tunneling 
coupling. We also build the BCS Hamiltonian in the band 
representation of the bilayer system; this is particularly 
useful to treat the interlayer hopping non-perturbatively. 
In Sec. IV, we thoroughly compute the superfluid stiff- 
ness and T c in the presence of a finite hole asymmetry. 




FIG. 1: Evolution of T c for the bilayer when both layers are 
characterized by prominent intrinsic d-wave correlations; Si 
and <$2 < 0.3 represent the hole densities of the two layers 
and (81 + 82) /2 is fixed to 0.2 close to optimal doping. The 
red dashed line is obtained from perturbation theory in the 
band basis. Parameters in the Hamiltonian (1) are J/t = 0.2, 
t' = and t±/t = 0.4. In this figure, T c , opt means T c for the 
bilayer system at optimal doping (see Fig. 5 for t±/t = 0.4). 



II. MODEL AND METHODOLOGY 

Our starting point is the renormalized low-energy the- 
ory [47l [62j which takes into account the proximity of the 
Mott insulating ground state. Essentially, the Gutzwillcr 
projector [63j ensuring that configurations with doubly 
occupied sites are forbidden is replaced by statistical 
weighting factors. The projection operator then is elimi- 
nated in favor of the reduction factor g t i = 25,-/(l + 5j) in 
the kinetic term [64j where Si represents the hole doping 
or the number of holes per site in the layer i — 1 or 2. 
In addition, the projection operator enhances spin-spin 
correlations in each layer: g S i — 4/(1 + Si) 2 [47l ]. 

The bilayer system in the strong interaction limit then 
is described by the general Hamiltonian: 



H = 



-t9tl £ 4o C 3o- ~ t'gtl £ ^ma C no + h.C 
<i,j>a <<m,n»a 

tg t 2 £ 4<T d i° ~ t'9t2 £ d\no d no + h.C. 
<i,j>cr 



<<m,n>>cr 



/'1 



£4 



M2 
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<i,j> <i,j> 

Jx 5s x£s?-Sf, (1) 
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where the operators c CT and d a represent electron opera- 
tors with spin a for layer 1 and 2, respectively, < i,j > 
and << m,n » refer to nearest-neighbor and next- 
nearest-neighbor pairs (and we have implicitly assumed 
i < j and similarly for m and n), and S c and S d de- 
note the spin-1/2 operators in each layer. In our model, 
the dopings of the two layers are independently tuned 
through the chemical potentials \i\ and pi. At a general 
level, the two layers are coupled via the single-particle 
tunneling contribution t±g t± and through the exchange 
term J±g s ±, where g t ± = ^Jgt\3t2 and g s± = ^/g s ig S 2- 

In Appendix A, we briefly introduce the methodology 
and the numerical procedure in the context of the sin- 
gle layer following Zhang and Rice [13] . The advantage 
of starting with this effective low-energy theory is that 
the d-wave superconducting ground state can be stud- 
ied essentially using the usual (unprojected) BCS wave- 
function. In the superconducting state, results obtained 
within this method are in excellent agreement with vari- 
ational Monte Carlo for projected d-wave states (48j . 

The main results for the single layer situation are 
presented in Fig. 2. We check that the Fermi liquid 

order parameter Xij — (3<7s>//4) Yla( c ia c j<r) * s almost 
doping-independent whereas the pairing order parame- 
ter A,j = (3g s J /4) y, , e aa > kc ia c, a i ) follows the pseu- 
dogap (crossover) line (the indices i and j here in- 
volve nearest neighbor sites); we look for mean- field so- 
lutions with Xij — Xj and Ajj = A along ^-direction and 
Aij = — A along y direction to ensure d-wave pairing. 
The d-wave nature of the order parameter here is dictated 
by t he p rominent antiferromagnetic fluctuations at (7r, tt) 
[23l |46| . Similar to Ref. [4!| here we assume a unique 
superconducting gap spreading over the full Fermi sur- 
face. In reality, the antinodal points of the Fermi surface 
are rather governed by the pseudogap [Tl|, In fact, 
we cannot exclude that the two-gap structure might arise 
from another competing order with the superconductiv- 
ity which may alter the results found below. On the other 
hand, the superfluid density can be formally derived from 
the quasiparticle contribution close to the nodal points. 
Hence, this argument rather supports that T c mostly de- 
pends on the superconducting gap. We have checked that 
our numerical approach to minimize the free energy per- 
fectly agrees with the mean- field equations (|A7I) . 

Within the "renormalized" mean-field theory (or 
equivalently the slave-boson theory [H, l65U67j p. the su- 
perfluid density at T = in Eq. (|B3I) is proportional to 
the hole doping 5, as confirmed experimentally (68j . 

The superconducting transition temperature of the 
isolated layer is evaluated using two complementary ap- 
proaches [43]. First, the renormalized mean- field theory 
predicts T c g t A (see Appendix B). The second ap- 
proach consists to evaluate the temperature dependence 
of the superfluid stiffness. A theory of T c for the under- 
doped cuprates has first been built by analogy to the 
Kosterlitz-Thouless transition in two dimensions. In- 
deed, Emery and Kivelson in 1995 proposed a model 
based on phase fluctuations [3] . In their picture, the pseu- 




FIG. 2: Magnitudes of the mean- field variables (in units of 
3g s J/4) and the zero-temperature superfluid density (in units 
of t) for the single layer described by the renormalized t—t' — J 
model versus hole doping. Here, we fix J/t = 0.2, t' = and 
T c is determined from T c ~ gtA. and from ps(T c ) = 0. 



dogap region is governed by phase fluctuations and at T c 
the superfluid stiffness jumps by the universal amount 
2T C / tt. On the other hand, as mentioned by Lee and Wen 
in 1997 [69] the thermal excitation of quasiparticles near 
the nodal points rather produce a linear decrease of the 
superfluid stiffness p s {T) with temperature. This is the 
earliest experimental evidence of d-wave symmetry [20| . 
Now, coming back to the Kosterlitz-Thouless scenario, 
this implies that the p s which controls the transition is 
not p s (0) but p s {T) which is greatly diminished by quasi- 
particle excitations, p s {T) ~ p s (0) — g{T) where g(T) is 
a linear function at low temperatures. Thus, T c can also 
be defined by p s (T c ) = 0. Then, T c can also be evaluated 
numerically using Eq. (Bl). Performing an expansion 
very close to the nodal points leads to [69j, [70( 



-> 0) = aT 

2 2 In 2 Vp 
1 , 

7T VA 



(2) 



and the ratio between the longitudinal and transverse 
velocities at the nodes reads (see Appendix B): 



v F /v A 



2tgt + x 
A ' 



(3) 



(We neglect the temperature dependence of x and A be- 
low T c as J ^> T.) The ratio vf/va is measured through 
the thermal conductivity [7l| . An important assumption 
made in Eq. ^ is that the d-wave quasiparticles are 
characterized by a renormalized current (70j — aevp (see 
Appendix B). We introduce the parameter a which can 
be seen as a phenomenological Landau parameter inher- 
ited from the normal state. In principle, the quasiparticle 
charge ae should be determined experimentally |49j . 
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To have a good agreement between the two definitions 
of T c and to reproduce the dome-shaped T C (S) phase dia- 
gram of the single layer (see Fig. 2), then we fix a ~ 0.9. 
Note that the value of a depends on the precise scheme 
used to treat interactions close to the Mott state [49(. 



III. BILAYER AT OPTIMAL DOPING 

First, we apply the methodology of Sec. II to the op- 
timally doped bilayer system. We intend to check that 
the superconducting state and therefore T c are rather 
stable toward single- (quasi)particle tunneling favored by 
the transverse hopping term tj_- In fact, since the trans- 
verse hopping term is still weakened by the Gutzwiller 
statistical weighting factor g t ± (which is equal to gt for 
symmetric layers), we shall show that T c is almost un- 
changed until moderate values of t± where gt±t±_ ~ J. 
Essentially, the prominent superconducting gap in each 
layer tends to prevent the proliferation of quasiparticles. 



A. Diagonalization in the Band Basis 

In the case of a symmetric bilayer model (with equal 
hole dopings Si = 5 2 = 5 and fii = fJ>2 = (J>) it is conve- 
nient to use the bonding and antibonding representation: 



bia /— {pier ~l~ d>io ) 

V2 



(4) 



which allows to diagonalize all the single-particle hopping 
terms and therefore to treat t± non-perturbatively. 

It is also judicious to introduce explicitly the mean- 
field order parameters for the two layers: 



Xij 


3 

<7 


Xij 


a 


Xii 






<T(T' 


A 2 


— 4*f? st ^ ^ ^ £<t<t' (dig dj(j! / 
00' 


A+ 





(5) 



For symmetric layers, we can define g S i = g s ± = g s = 
4/ (1 + S) 2 , and we look for mean-field solutions xjj = Xij 

Xij = X2 where \i = X2 = X, Xu = XX, A^ = +A X for 
two nearest neighbors along x direction and — Ai for two 



nearest neighbors along y direction, and similarly for the 
second layer with Aj = A 2 = A. 

We also check that for < t±/t < 0.5, the order pa- 
rameter A^ is always negligible which means that the 
only pairing contribution is the intralayer pure d-wave 
contribution. (In this paper, we are not interested in 
the regime of (very) large interlayer transverse hopping 
amplitudes.) Hereafter, we thus omit the negligible con- 
tribution from A^. Therefore the main coupling be- 
tween the layers is the single-(quasi)particle tunneling 
term (when assuming J± cx t\/U <C t±) and the term 
J± just renormalizes t± by producing a finite \±- 

In the band basis, the mean-field Hamiltonian reads: 



Hxin + Ha + H, 



const ; 



(6) 



where 



H 



Hxin — E $k,&frk CT frkq + E £k,qQk ff ^kg 
ko- 

k 

- N. E 



-14 



aU, kt u -H 



h.c. 



i=l,2,± 



I A, 



2N S ^5. (7) 



Our convention for the chemical potential follows that of 
Ref. [53] and N s is the total number of sites. We identify: 

£,k,b/a = -(2tg t + x)[cos(k x ) + cos{k y )} 

- [4g t t' cos(k x ) cos(k y ) ±g t t± ± x±] - M- ( 8 ) 

For simplicity, the lattice spacing is set to unity and for 
symmetric layers, g t ± = VStiS^ = 9t = 26/ (1 + S). The 
Fermi surfaces associated with the two bands get splitted 
as a result of the transverse hopping amplitude t± and 
X±- Further, the pairing parameters Ak, a and Ak,b are 
coupled through the mean-field order parameter A; more 
precisely, neglecting A^ results in 



A k , Q = A k ,6 = A(cos(k x ) - cos(fey)). 



(9) 



Interestingly, one can easily diagonalize H sym for any 
value of tj_ and the mean-field free energy is given by 
(again xi = X2 = X, &i = A 2 = A and A ± = 0): 



F. 



MF 



-2T ^ 



cosh 



N, 



E 

i=l,2,_ 



\Xi 



Ek t i 
~2T 

A,; 1 2 



2N sf j,6. (10) 



The quasi-particle excitation energy for each band is rem- 
iniscent of the BCS theory: 



Ek, a /b — V (£k,a 



7,/b) 



(A 



k,a/b) 



(11) 



The mean-field equations can be obtained by minimizing 
the free energy with respect to /i, Xi A and x± along the 
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FIG. 3: Magnitudes of the order parameters (in units of 
3g s J/4) and superfiuid density (in units of t) for the opti- 
mally doped bilayer system at zero temperature as a function 
of t±. Optimal doping means S = 0.17, J/t — 0.2 and for this 
figure, t'/t = -0.3 / 0. 



lines of the single layer case and at zero temperature: 



1 \ \ ^ 6m 

2N a A 4- E k , 



(12) 



i—a,b k 



^~9sJ Y(cos(k x ) + cos(k y )) 



,.b k 



2E k . 



A = 



X-L 



87V. 



2-Ek ,i 



! 5 s J X! E( cos ( fc:c ) ~~ cos(fc y )) Ak 

i—a : b k 
j—a,b k J 



where in the last line j — + when j = a and j = — when 
j = b. Results for the mean-field order parameters at 
zero temperature are shown in Fig. 3. Here, it is worth 
mentioning that even though x± increases as a result 
of the finite ij_ the d-wave gap remains almost constant 
reflecting the stability of the d-wave state toward the 
intcrlayer single-(quasi)particle tunneling term. 

Further, it should be noted that in the band represen- 
tation the two bands are still characterized by the same 
d-wave gap. However, details of their distinct Fermi sur- 
faces may matter when evaluating the superfiuid stiffness. 



B. Superfiuid density and T c 

The diamagnetic current is {Jd)^i = D^ V A V where 



ktr 



dkndk v 



(4a C k<7 + 4<Aa 



(13) 



and is the "bare" spectrum (the magnetic term J does 
not contribute to the electric current): 

£k = —2tg t (cos(k x ) + cos(ky)) — 4g t t' cos(k x ) cos(k y ) . 

(14) 

The vector potential A is directed along the layers. Then, 
we can rewrite D^ v in the band basis as: 



ka 



dkndk v 



(15) 



Using the standard BCS ground state wavefunction and 
the mean-field Hamiltonian in the band representation, 
then we identify: 



ka) = 1 — — — tanh 



Ek.b 



V 2T 



(16) 



We obtain (a^a^) in a similar way. The zero tempera- 
ture part of the superfiuid density then takes the form: 



k,i— a,b 



£k,i 
Ek,, 



(17) 



Here, ^ = d 2 (Q/dk 2 means that the vector potential 
is directed along the x-axis. For small t±, we find the 
following expression: 



Ps(T = 0) 
2p s (0) 



1 



2 P s(0) 



k 



2 ^ k ^k.O 



(18) 



with p s (0) being the zero-temperature superfiuid den- 
sity of the isolated layer and £k an d E k of the isolated 
layer are defined in Appendix A. At a general level, 
Sk£k.o£k/^k < which tends to suggest a light down- 
turn of the zero-temperature superfiuid density when 
switching on the interlayer coupling t±. On the other 
hand, we check that the exact superfiuid density at zero 
temperature computed from Eq. (Bl) does not substan- 
tially decrease until quite large values of t± ~ 0.5t. At 
small values of t±, we note an excellent agreement be- 
tween the exact expression of p s (T = 0) in Eq. (Bl) and 
the weak-coupling approximation in Eq. (|18[) ; see Fig. 4. 

To compute the critical temperature T c for the bilayer 
system at optimal doping, first we follow Goren and Alt- 
man (49j and diagonalize the mean-field Hamiltonian in 
the layer basis. The critical temperature T c is defined 
via Eq. (Bl) by p s (T c ) = in the layer and band ba- 
sis. Remember that the transition temperature T c can 
be equivalcntly defined in the band basis using Eq. (10). 

At low temperatures, the linear-T dependence of the 
superfiuid stiffness essentially stems from the paramag- 
netic component (see Appendix B). In the band basis, 
interestingly, this can be separated into bonding and anti- 
bonding contributions. Close to the nodal points, we get 



p 5 (T)« Ps (T = 0) 



i—a,b 



21n2 

7T V A , 



(19) 



6 



P s (T=0) 



1.05 



□ — □ analytical approx 
3fc — & numerical 



Tc / T 




FIG. 4: p a (T = 0) from Eq. (Bl) (in units of t) and a weak- 
coupling approximation in the band basis (t 1 — 0). 



where formally 



FIG. 5: Transition temperature for the bilayer at optimal dop- 
ing, normalized to T c of the isolated layer for the same doping 
(denoted T c o in this figure), versus t±. The two curves corre- 
spond to the two approaches defined in Sec. Ill for accessing 
T c . We use the parameters t' = 0, J/t = 0.2 and 8 = 0.17. 



Vf 
t 

VFA 



(20) 



IV. ASYMMETRIC LAYERS 



and vp for the single layer has been defined in Appendix 
B. It is relevant to mention that in the band basis close 
to the nodal points the longitudinal and transverse ve- 
locities obey: 



VF,a v F ,b 2tg t + X 



VA,a 



VA.b 



(21) 



This shows that the ratio {vF, a /b/ v A,a/b) remains quite 
constant until moderate values of t±; in particular, it 
does not involve x ± ■ This allows us to safely approximate 
on ~ a ~ 0.9. We shall also mention that since the ratio 
VF,i/vA.i remains almost identical to that of the isolated 
layer until moderate values of t± this already suggests a 
very slow reduction of T c with t± . 

For small values of t±, the zero-temperature value of 
the gap obeys: 



A(tj 



(t±9t + X±) 
4£3 



(cos(fc x ) — cos(k y )Y 



(22) 

(where A corresponds to the value of the gap for the 
isolated layer at optimal doping; see Fig. 2). For small 
tj_) we predict T c — T c q — C(gtt± + x±) 2 with C > 0. 

The curve of T c versus t± for the optimally doped bi- 
layer system is presented in Fig. 5. This unambiguously 
confirms that the bilayer system at optimal doping is 
rather stable toward single-(quasi)particle tunneling un- 
til quite large values of t±, defined roughly by g t t± ~ J 
(J controls the pairing properties of the two layers). The 
prominent superconducting gap hinders the proliferation 
of quasiparticles close to the nodal points. 



In this Section, we address the case of asymmetrically 
doped layers, i.e., <C 5\ < 0.2 and 0.2 < S 2 < 0.3 
(0.2 is roughly the hole density in each layer at optimal 
doping). We seek to understand if such a finite asym- 
metry in the hole dopings of the layers will result in an 
increase or in a decrease of T c . Using the band basis, we 
show how a finite hole asymmetry will result in a pairing 
term coupling the bonding and antibonding bands which 
helps diminish the quasiparticle nodal contribution for 
weak asymmetries, then enhancing T c of the optimally 
doped situation. The results derived below assume that 
each layer exhibits intrinsic d-wave pairing correlations. 



A. Interband pairing term 



More precisely, the Hamiltonian becomes H = H syr , 
H asy where H sym can be found in Sec. Ill and 



H aS y 



kcr 



E A k (4t&-ki + & kt«-k, + h.c). (23) 



k 



In this case, the c-electrons of the first layer and the d- 
electrons of the second layer are characterized by different 
band structures £k,i and £k,2 respectively (since gn 7^ gt2 
and g s i ^ g 8 2)- In general, one can decompose 



6c,i 

£k,2 



fat) I sas 
?k + ?k 



fai? fas 
Sk — <5k ■ 



(24) 
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Then, we identify: 



(25) 



*(ffti - gti) + ^ (Xi - X2) ) (cos^) + cos(fc y )) . 



For simplicity, hereafter we assume that t' = 0. Similarly, 
one can rewrite the pairing order parameters of the two 
layers as 



A k ,i = A^ + A£ s 
A k , 2 = A^-A£ s , 

which also results in 



(26) 



A£ s = -(A 1 -A 2 )(cos(fc x )-cos(fc !/ )). (27) 



Concerning the H sym part this involves and A£": 



6k,6/a 



(28) 



%ti + .9*2) + ^(Xi + X2) ) (cos(fc x ) + cos(fc y )) 



A k , a = A k:6 = -(Ai + A 2 )(cos(k x ) - cos(ky)). 

One can still compute the critical temperature T c and 
the superfluid density p a {T c ) (from Eq. (Bl)) in the layer 
basis by diagonalizing the Hamiltonian. On the other 
hand, to gain some intuition, we also treat the asymme- 
try terms in the Hamiltonian to second order in perturba- 
tion theory (which is justified for quite small asymmetries 
around the optimal doping). 

Essentially, in the band basis, this results in correc- 
tions to £k,a/& and A k , a /6 such that £k,b and A k; b become 
£ k> b and A k ,fc defined as (see Appendix C): 

- ((er) 2 Ck, Q - (An 2 6c, a + ^ k s A k , a ) (29) 



A 



k,6 



v k,a 
1 



((A-) 2 A k , a - (?r) 2 A k ,a + A^k.a) , 



and similarly for £ k , a and A ki<l . (There is no first order 
correction to the ground state energy.) Note that even 
though the main pairing contribution is the pure d-wave 
intraband component, the interband pairing correlations 
favored by A£ s will contribute to reduce the nodal quasi- 
particle contribution to the superfluid stiffness p s {T). 

From second-order perturbation theory, the effect of 
££ s and A£ s is primarily to renormalize the band struc- 
ture parameters entering into the energies E^ a and _E kj b 
of the quasiparticles at the nodal points. 



B. Renormalization of nodal contributions 

The next step to compute the transition temperature 
T c is to check how £ kja /f, and A k a /f, vary to linear order 
in k-space for points around the nodal points. 



Notice that the bonding (b) Fermi "surface" is larger 
than the antibonding (a) Fermi surface. This is a very 
general fact stemming from the finite interlayer hopping 
term t±. Therefore, we can denote Mb and M a the nodes 
associated with the bonding and antibonding bands re- 
spectively and k a f, > (proportional to t±) is the k-space 
distance between the bonding and antibonding Fermi sur- 
faces. Now, (fci,fc 2 ) are local coordinates with origin at 
Mb such that: 

£k,6 = v F ,bh (30) 
A k ,6 = v A ^ b k 2 

£k,a = « F ,«(fel + k ab ) 

A k:Q = v A , a k 2 , 

where ki,k 2 — > 0. Then, terms such as (A£ s ) 2 are of 
order k 2 and therefore do not contribute to linear order. 
As a result, we can approximate 

(er) 2 £k, Q 



£k,6 W £k,b - 



El 
k,a 



(31) 



i 



k,6 



VF,akab 



Kab 



The first correction term in £ ki b shifts the position of the 
node Mb and the second term gives a correction to VF,b- 



VF,b = v F h ■ 



f 



VF,ak 2 



(32) 



a 6 



Writing A£ s = v A s k 2 we also get: 



1 



^k,fc 



^k,fc- 



{vF.akab)'- 



Therefore, v A ,b is renormalizcd as: 
1 



(2^ B V F .akabV A S ~ (SvU Va) fe- 

(33) 



V A ,b ~ «A,6- 



{vF.akab) 2 

Similarly, we also find 



{2CP B V F ,ak ab V A S 



and, 



V A ,a ~ VA,a- 



V F ,a = V F ,a 



1 



v F,bk 2 ab 



) V«) • 
(34) 



(35) 



{vF,bkab) 



(36) 

Now, a close inspection of all the contributions leads to: 



VF,a 



V F ,b 



V A M v At b 



+ 



VF,a 
V A ,a 
VF,b 



A 



1 + 



V A , a VF,bkal 
as n .as 
b V A 



(37) 



V A ,b \ VA,bVF,akab, 

It should be noted that for very small values of t± such 
that k a b — > then we could approximate £™ ~ £™ 
and the nodal corrections would have practically no ef- 
fect. For the symmetric case, remember that vp.a/v Aya = 
VF,b/v A ,b and we also check that v Aia VF.b = v Ai bVF.a- 
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C. Enhancement of T c 

Now, let assume a moderate value of t± such that 
<C gt±t± -C J- At a general level, we can compute 
£jy and £jy numerically and then extract the supercon- 
ducting transition T c . Here, it should be noted that since 
a prominent tj_ makes cos(k x ) + cos(k y ) larger for a small 
Fermi surface and (gti — g-ei) < then this immediately 
implies £%f A > £ff B (for /xi ~ [i 2 and xi ~ Xa)< Based on 
the nodal quasiparticle contribution only, in the case of 
a finite (small) hole asymmetry, one then predicts: 

^«(l + ^A^A W^A \ >1(38) 

Tc.opt V VA,bVF,akab VA, a VF,bk a b J 

Here, T CjOJrt refers to T c at optimal doping for the sym- 
metric bilayer. In Fig. 1, we compare the result of Eq. 
(I38p valid for a small doping asymmetry with the nu- 
merical results obtained in the layer basis directly. Both 
results seem in excellent agreement for small asymme- 
tries. This traduces that additional (weak) supercon- 
ducting fluctuations mediated by A£ s mostly affect the 
nodal contribution(s) whereas other regions (in fc-space) 
are protected by the large d-wave gap. 

The enhancement of T c is attributed to the interplay 
between interband pairing correlations {v A B ) and the fi- 
nite interlayer hopping. This result relies on the exis- 
tence of a large and small Fermi surface induced by a 
finite interlayer hopping term. In Fig. 1 we maintain the 
average hole density ((5i+<5 2 )/2 ~ 0.2 fixed and each layer 
exhibits intrinsic d-wave pairing correlations (the param- 
eter a ~ 0.9 for each layer). In Fig. 6, in contrast, we 
consider an overdoped layer described by a free-electron 
like model (with a.2 = 1) and essentially we corroborate 
the result obtained in Ref. [49], choosing exceptionally 
a.\ = 0.5 for the underdoped layer. Hence, we conclude 
that two scenarios allow to increase T c in bilayer sys- 
tems, one relies on the presence of a very overdoped free- 
electron layer and the other relies on the enhancement of 
pairing fluctuations in the band basis induced by a finite 
hole asymmetry around optimal doping. 

V. CONCLUSION 

In this paper we corroborate that significant enhance- 
ment of T c in strongly-correlated heterostructures is pos- 
sible under realistic conditions. More precisely, apply- 
ing an effective low-energy theory of the doped Mott in- 
sulator we have investigated a bilayer Hubbard system 
where both layers exhibit intrinsic intralayer (intraband) 
d-wave superconducting correlations. Using the renor- 
malized mean- field theory which is usually well controlled 
when focusing primarily on the superconducting state, 
we have evidenced that the increase of T c results from 
the delicate balance between the moderate single-particle 
tunneling term coupling the layers and the finite hole 
asymmetry around optimal doping which tends to reduce 



°- 04 i 




FIG. 6: Here, we consider that the overdoped layer is suffi- 
ciently doped such that it is described by a free-electron model 
(a.2 = 1 and Xi ~ 0) 5H> * ne doping level of the metallic layer 
is &2 = 0.35. Parameters are t±/t = 0.5, J/t = 0.3 and similar 
to Ref. [49|, for this figure, exceptionally we set the renor- 
malized quasiparticle charge of ai = 0.5 and <5i = 5. Here, 
T^ ax refers to the maximal value of T c for the single layer. 
The enhancement of T c is in agreement with Ref. [49| (but it 
clearly depends on the precise choice of ati.) 

the quasiparticle contribution to the superfluid stiffness 
by reinforcing superconducting fluctuations. In fact, we 
have built the BCS Hamiltonian in the band representa- 
tion of the bilayer system which is particularly judicious 
to build a non-perturbative theory in the single-particle 
tunneling term coupling the layers. We have also shown 
that the d-wave superconducting state is quite robust to- 
ward the interlayer single-(quasi)particle tunneling term. 
The key point to enhance T c in these heterostructures 
is that a finite interlayer hopping produces a larger and 
smaller Fermi surface and a moderate hole asymmetry 
between the two layers reinforces (interband) supercon- 
ducting fluctuations. This scenario requires that both 
bands are filled (vp, a =^ and vp^ ^ 0) and hence t± 
should not be too large; a too large t± rather favors single 
quasi-particle interlayer tunneling and therefore is gener- 
ally not helpful for superconductivity (see Fig. 5). 

It is important to distinguish this scenario based on 
two layers exhibiting prominent d-wave pairing from the 
other scenario based on a very overdoped layer described 
by a metallic bath which essentially serves to increase 
the number of careers [4£| (see also Fig. 6). The latter 
situation using a very overdoped metallic layer seems to 
lead to a more substantial increase of T c (note that the 
values of a\ in Figs. I and 6 are different). 

It should also be mentioned that in this paper we have 
ignored gauge (phase) fluctuation effects which may play 
an important role for (quasi-) two-dimensional systems. 

We thank T. M. Rice, C. Ahn and F. Walker for 
stimulating dicussions. This work is supported by 
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Appendix A: Doped Mott insulator and d-wave 
superconductivity 

The effective Hamiltonian of a doped-Mott insulator 
takes the form: 



Ht-t'-j = -tg. 

- t g t 



+ h.c. 



^ t c \na c na 
<<m,n>>er 



<i,3> 



h.c. 



S?-S5. (Al) 



Here, < i,j > represent the nearest-neighbor pairs and 
<< m, n » the next nearest-neighbor pairs (and we 
implicitly assume that i < j and similarly for m and n) 

lg and spin 

LV > ] TT^m 

respectively and 6 is the hole doping. Below, we closely 
follow the notations of Zhang and Rice [47| . 

We then introduce the mean-field order parameters: 



The statistical weighting factors for the hopping and spin 
exchange coupling are g t = ^ ^64] and g s = jj^p 



The mean-field free energy is given by 

in cosn i 



Fmf 



-2T$> 



2T 



\\\ 



IAI 



\g s J \g s J 



N s /j,S (A5) 



where the quasi-particle excitation energy obeys 



(A6) 



The Boltzmann constant is set to unity. 

The mean-field equations can be obtained by directly 
minimizing the free energy with respect to x, A and by 
imposing (OFmf /d/i) = 0. The mean-field variables of 
Fig. 2 are solutions of the following equations: 



5 = — tanh ( — ^ 

N s V V 2T 

k 



x = 



4AT, 



s J^2(cos(k x ) + cos tanh 



2T 



k ^ \ 



(A7) 



Appendix B: Superfluid Density and T c 



Xij 



2 ^ 



<J(T ! 



(A2) 



with e-^i = 1 = — £4.^ and zero otherwise. Then, we look 
for mean-field solutions with Xij = Xj an d A,-j — A along 



x-direction and A.; 



-A along y direction to ensure d- 



wave pairing. The mean-field Hamiltonian then reads: 



H, 



MF 
t-t'-J 

HRin 

Ha 

Hronst 



HKin + Ha 



H r , 



krr 



Al 2 



h.c. 



k 



+ 



l<? s J ' §<? S J 



(A3) 



where A s is the total number of sites (and the chemical 
potential has been introduced following Zhang and Rice 
|47|). Assuming a square lattice geometry, we identify: 



£ k = -(2tg t + x)[cos(fe x ) + cos(ky)] 

— Ag t t' coa(k x ) cos(k y ) — fj, 
Ak = A(cos(fc x ) — cos(kyj). 



(A4) 



The superfluid density is formally defined as: 



1 
Vol 



d 2 F, 



MF 



dAf,dA v 



(Bl) 



A=0 



where Vol is the volume of the system and A is the vector 
potential entering in the kinetic part of the Hamiltonian 
through a phase: 



JeA; 



cLc-, 



H K in(A) = -tg t Y 

<i,j>0 

t 9t ^ e ***" c lna c no 



<<m,n>><7 



(B2) 



Here, without loss of generality we assume the vector po- 
tential to be along the a;— axis: A = A x . The superfluid 
density at any temperature can therefore be evaluated 
numerically through Eq. (|F31|) . 

Further, at T = 0, the superfluid density can be easily 
obtained analytically, and it is given by (see Sec. Ill B): 



(o) = ^E&',o(i-&/£k), 



(B3) 



where f£ = d 2 Q/dk 2 x with £° = -2tg t (cos(k x ) + 
cos(ky)) — 4g t t' cos(k x ) cos(fc y ) being the hopping part of 



the kinetic energy. [Below, we neglect the T-dependence 
of this diamagnetic contribution since a power-counting 
argument shows that this T-dependence is 0(T 2 ).] 

At finite temperature, the superfluid density is in- 
evitably suppressed by the normal state quasiparticle ex- 
citations near the four nodal points q = (±g, ±q) with 
q = ir/2 at half-filling. In the vicinity of the node (q,q), 
we have the anisotropic Dirac spectrum: 



E k 



v F k\ 



u A n 2 



(B4) 



where for the square lattice 



Vp 
Va 
cos(g) 



V2(2tg t + X ) sin(<?) 
\/2A sin(g) 

2(2tg t + x) ' 



(B5) 



vf and va are the nodal quasiparticle velocities in the 
longitudinal and transverse directions, respectively. For 
simplicity, here we assume that t' = allowing a simple 
analytical solution. More precisely, by definition 



v A = 



Vp 



dk 2 
dE k 
dki 



|fci,fc 2 ->0 



|fcl,fc 2 ->0 



(B6) 



where fci = (k x + k y — 2q)/\/2 and k 2 = (k x — ky)j\f2 with 
q being associated with the Fermi momentum kp = (<y, q) 
on the nodal point (q — it/2 at half- filling on the square 
lattice if t' = 0). The Fermi momentum obeys: 

fk, = = -{2tg t + X )(cos(k x ) + cos(k y )) - p. (B7) 



Therefore, we get cos(g) = 
and Ak near a given node, 



2(2t gt +x) ■ 



We can expand £k 



£ k w V2(2tg t + x)sm(q)h 
A k « -V2Asin(g)fc 2 - 



(B8) 



This results in Eqs. (jB4j) and p5|) . 

In the presence of a vector potential, the quasiparticle 
spectrum exhibits a shift: 



£(A,A) = £(fc)-j(fc)-A, 



(B9) 



where the current j carried by the normal state quasi- 
particles can be formally written as: 



-eavp. 



(BIO) 



Ignoring interactions between the quasiparticles would 
result in a = gtvp/vp where v F corresponds to the bare 
Fermi velocity when setting x — 0- Such a choice of 
a would not allow to reproduce the dome-shaped T c (§). 
Therefore, in this paper, a ~ 0.9 will be rather taken as 
an effective (constant, doping independent) parameter 
[70j which can also be regarded as an effective charge. In 




0— © numerics 
— analytical approximation 



FIG. 7: ps{T) in units of t using Eq. (Bl). The analyt- 
ical approximation (Bll) corresponds to a low-temperature 
expansion very close to the nodes. The parameters are the 
same as in Fig. 2 and a is fixed around 0.9. Here, 5 = 0.17. 



Eq. (Bl), we adjust the vector potential A such that it 
reproduces the correct effective charge. 

Using Eqs. dA"5|) . pi). ([B4| and (|B9]l. and perform- 
ing the integral in Eq. (|A5j) in momentum space near the 
nodes gives the low-temperature linear approximation: 



Ps(T) 



Ps(T = 0)- 
2 2 In 2 Vp 

7T VA 



aT 



a 



(Bll) 



(B12) 



The ratio va/vf near a nodal point reads: 
Vf _ 2tg t + x 
v A ~ A 

The transition temperature T c is determined through 
Eq. (Bl) as the temperature at which the superfluid 
density vanishes due to the thermal excitation of quasi- 
particles (but, not necessarily nodal); see Fig. 7. 

Setting a ~ 0.9 in Eq. (Bl) through the vector po- 
tential A allows us to recover a form of T c {8) which is 
reminiscent of the superconducting order parameter 



A 2 
^sc 



(BCS\P(4 t clc l+n c ]+li )P\BCS) ~ 5t 2 A 2 . 

(B13) 

for a large distance I, as shown in Fig. 2. 

Appendix C: Perturbation Theory 

Here, we derive Eqs. (29) in the main text obtained 
by treating H asy in perturbation theory. First, it is con- 
venient to write the Hamiltonians of each band as 2 x 2 
block matrices. When H asym = then this results in 



Ha 



6 
A 6 

A a 



A ft 
A a 



(CI) 
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For simplicity, we suppress the momentum index k. The 
mixing between the bonding and antibonding sectors is 
given by V AB = V BA = H asy : 



Hasy — 



£>as 

A a , 



(C2) 



Solving the time-independent Schrodinger equation and 
integrating out the A-subsystem the effective Hamilto- 
nian for B subsystem is given by: 



Hb = H B + Vba(-Ha) 1 Vab 
with V B a{-H a )- 1 Vab = 



1 



£,as 



-A„. 



-A a 



A„ 



A as 



(C3) 



(C4) 



Therefore, we check that 



Hb + Vba(-Ha) 1 Vab = 



6 
A„ 



A b 

-lb 



(C5) 



where and A{, are precisely defined as 



6 = 6 - ^ ((Z aS ) 2 ta - (A QS ) 2 £ Q + 2r s A as A a )(C6) 



1 

m 



A b = A b - ((A as ) 2 A a - (r s ) 2 A a + 2£ as A as Q 
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